Targeting G9a translational mechanism of SARS-CoV-2 pathogenesis for multifaceted therapeutics of COVID-19 and its sequalae (preprint)

By largely unknown mechanism(s), SARS-CoV-2 hijacks the host translation apparatus to promote COVID-19 pathogenesis. We report that the histone methyltransferase G9a noncanonically regulates viral hijacking of the translation machinery to bring about COVID-19 symptoms of hyperinflammation, lymphopenia, and blood coagulation. Chemoproteomic analysis of COVID-19 patient peripheral mononuclear blood cells (PBMC) identified enhanced interactions between SARS-CoV-2-upregulated G9a and distinct translation regulators, particularly the N6-methyladenosine (m6A) RNA methylase METTL3. These interactions with translation regulators implicated G9a in translational regulation of COVID-19. Inhibition of G9a activity suppressed SARS-CoV-2 replication in human alveolar epithelial cells. Accordingly, multi-omics analysis of the same alveolar cells identified SARS-CoV-2-induced changes at the transcriptional, m6A-epitranscriptional, translational, and post-translational (phosphorylation or secretion) levels that were reversed by inhibitor treatment. As suggested by the aforesaid chemoproteomic analysis, these multi-omics-correlated changes revealed a G9a-regulated translational mechanism of COVID-19 pathogenesis in which G9a directs translation of viral and host proteins associated with SARS-CoV-2 replication and with dysregulation of host response. Comparison of proteomic analyses of G9a inhibitor-treated, SARS-CoV-2 infected cells, or ex vivo culture of patient PBMCs, with COVID-19 patient data revealed that G9a inhibition reversed the patient proteomic landscape that correlated with COVID-19 pathology/symptoms. These data also indicated that the G9a-regulated, inhibitor-reversed, translational mechanism outperformed G9a-transcriptional suppression to ultimately determine COVID-19 pathogenesis and to define the inhibitor action, from which biomarkers of serve symptom vulnerability were mechanistically derived. This cell line-to-patient conservation of G9a-translated, COVID-19 proteome suggests that G9a inhibitors can be used to treat patients with COVID-19, particularly patients with long-lasting COVID-19 sequelae.

A Murine Model of Post-acute Neurological Sequelae Following SARS-CoV-2 Variant Infection (preprint)

Viral variant is one known risk factor associated with post-acute sequelae of COVID-19 (PASC), yet the pathogenesis is largely unknown. Here, we studied SARS-CoV-2 Delta variant-induced PASC in K18-hACE2 mice. The virus replicated productively, induced robust inflammatory responses in lung and brain tissues, and caused weight loss and mortality during the acute infection. Longitudinal behavior studies in surviving mice up to 4 months post-acute infection revealed persistent abnormalities in neuropsychiatric state and motor behaviors, while reflex and sensory functions recovered over time. Surviving mice showed no detectable viral RNA in the brain and minimal neuroinflammation post-acute infection. Transcriptome analysis revealed persistent activation of immune pathways, including humoral responses, complement, and phagocytosis, and reduced levels of genes associated with ataxia telangiectasia, impaired cognitive function and memory recall, and neuronal dysfunction and degeneration. Furthermore, surviving mice maintained potent T helper 1 prone cellular immune responses and high neutralizing antibodies against Delta and Omicron variants in the periphery for months post-acute infection. Overall, infection in K18-hACE2 mice recapitulates the persistent clinical symptoms reported in long COVID patients and may be useful for future assessment of the efficacy of vaccines and therapeutics against SARS-CoV-2 variants.

Resistance mechanisms of SARS-CoV-2 3CLpro to the non-covalent inhibitor WU-04 (preprint)

Drug resistance poses a significant challenge in the development of effective therapies against SARS-CoV-2. Here, we identified two double mutations, M49K/M165V and M49K/S301P, in the 3C-like protease (3CLpro) that confer resistance to a novel non-covalent inhibitor, WU-04. Crystallographic analysis indicates that the M49K mutation destabilizes the WU-04 binding pocket, impacting the binding of WU-04 more significantly than the binding of 3CLpro substrates. The M165V mutation directly interferes with WU-04 binding. The S301P mutation, which is far from the WU-04 binding pocket, indirectly affects WU-04 binding by restricting the rotation of 3CLpros C-terminal tail and impeding 3CLpro dimerization. We further explored 3CLpro mutations that confer resistance to two clinically used inhibitors: ensitrelvir and nirmatrelvir, and revealed a trade-off between the catalytic activity, thermostability, and drug resistance of 3CLpro. We found that mutations at the same residue (M49) can have distinct effects on the 3CLpro inhibitors, highlighting the importance of developing multiple antiviral agents with different skeletons for fighting SARS-CoV-2. These findings enhance our understanding of SARS-CoV-2 resistance mechanisms and inform the development of effective therapeutics.

Less neutralization evasion of SARS-CoV-2 BA.2.86 than XBB sublineages and CH.1.1 (preprint)

The highly mutated BA.2.86, with over 30 spike protein mutations in comparison to Omicron BA.2 and XBB.1.5 variants, has raised concerns about its potential to evade COVID-19 vaccination or prior SARS-CoV-2 infection-elicited immunity. In this study, we employ a live SARS-CoV-2 neutralization assay to compare the neutralization evasion ability of BA.2.86 with other emerged SARS-CoV-2 subvariants, including BA.2-derived CH.1.1, Delta-Omicron recombinant XBC.1.6, and XBB descendants XBB.1.5, XBB.1.16, XBB.2.3, EG.5.1 and FL.1.5.1. Our results show that BA.2.86 is less neutralization evasive than XBB sublineages. Among all the tested variants, CH.1.1 exhibits the greatest neutralization evasion. In comparison to XBB.1.5, the more recent XBB descendants, particularly EG.5.1 and FL.1.5.1, display increased resistance to neutralization induced by parental COVID-19 mRNA vaccine and a BA.5-Bivalent-booster. In contrast, XBC.1.6 shows a slight reduction but remains comparable sensitivity to neutralization when compared to BA.5. Furthermore, a recent XBB.1.5-breakthrough infection significantly enhances the breadth and potency of cross-neutralization. These findings reinforce the expectation that the upcoming XBB.1.5 mRNA vaccine would likely boost the neutralization of currently circulating variants, while also underscoring the critical importance of ongoing surveillance to monitor the evolution and immune evasion potential of SARS-CoV-2 variants.

A single-dose of intranasal vaccination with a live-attenuated SARS-CoV-2 vaccine candidate promotes protective mucosal and systemic immunity (preprint)

An attenuated SARS-CoV-2 virus with modified viral transcriptional regulatory sequences and deletion of open-reading frames 3, 6, 7 and 8 ({triangleup}3678) was previously reported to protect hamsters from SARS-CoV-2 infection and transmission. Here we report that a single-dose intranasal vaccination of {triangleup}3678 protects K18-hACE2 mice from wild-type or variant SARS-CoV-2 challenge. Compared with wild-type virus infection, the {triangleup}3678 vaccination induces equivalent or higher levels of lung and systemic T cell, B cell, IgA, and IgG responses. The results suggest {triangleup}3678 as an attractive mucosal vaccine candidate to boost pulmonary immunity against SARS-CoV-2.

S:D614G and S:H655Y are gateway mutations that act epistatically to promote SARS-CoV-2 variant fitness (preprint)

SARS-CoV-2 variants bearing complex combinations of mutations have been associated with increased transmissibility, COVID-19 severity, and immune escape. S:D614G may have facilitated emergence of such variants since they appeared after S:D614G had gone to fixation. To test this hypothesis, Spike sequences from an immunocompromised individual with prolonged infection, and from the major SARS-CoV-2 variants of concern, were reverted to the ancestral S:D614. In all cases, infectivity of the revertants was compromised. Rare SARS-CoV-2 lineages that lack S:D614G were identified and the infectivity of these was dependent upon S:Q613H or S:H655Y. Notably, Gamma and Omicron variants possess both S:D614G and S:H655Y, each of which contributed to infectivity of these variants. All three mutations, S:Q613H, S:D614G, and S:H655Y, stabilized Spike on virions, consistent with selection of these mutations by a common molecular mechanism. Among sarbecoviruses, S:Q613H, S:D614G, and S:H655Y are only detected in SARS-CoV-2, which uniquely possesses a polybasic S1/S2 cleavage site. Results of genetic and biochemical experiments here demonstrated that S:D614G and S:H655Y are likely adaptations to the cleavage site. CryoEM revealed that both mutations shift the Spike receptor binding domain towards the open conformation required for ACE2-binding and Spikes bearing either S:D614G or S:H655Y spontaneously mimic the smFRET signal that ACE2 induces in the parental molecule. Data from these orthogonal experiments demonstrate that S:D614G and S:H655Y are convergent adaptations to the polybasic S1/S2 cleavage site, which stabilize S1 on the virion in the open RBD conformation that is on-pathway for target cell fusion, and thereby act epistatically to promote the fitness of variants bearing complex combinations of clinically significant mutations.

HIF-dependent induction of alveolar miR-147 dampens SARS-CoV-2 immune evasion (preprint)

Acute respiratory distress syndrome (ARDS) causes significant morbidity and mortality during severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections1,2. Nevertheless, most patients with coronavirus disease 2019 (COVID-19) recover seamlessly without developing ARDS3, suggesting the existence of endogenous pathways to protect the lungs. Since many microRNAs (miRNAs) serve important roles in endogenous regulatory pathways, we pursued functional roles of miRNAs in COVID-19-ARDS. A screen of miRNAs in human alveolar epithelia or mice infected with SARS-CoV-2 identified miR-147 as a lead candidate. Transcriptional analysis implicated hypoxia-inducible factor 1A (HIF1A) in miR-147 induction during alveolar injury or SARS-CoV-2 infection. Functionally, mice with alveolar epithelial deletion of miR-147 showed increased lung injury in response to SARS-CoV-2 infection. mRNA sequencing and subsequent in silico miR-147-target analysis revealed reduced antiviral responses and identified SARS-CoV-2 ORF8 as a direct miR-147 target. Moreover, mice infected with SARS-CoV-2 and treated with miR-147 packaged in DOPC nanoliposomes showed significant protection with enhanced antiviral responses and improved survival. Finally, proof-of-principle studies in patients with COVID-19 highlighted this pathway in human SARS-CoV-2-associated ARDS. Together, our findings identify a previously unrecognized role of miR-147 in lung protection during viral pneumonia by attenuating immune evasion of SARS-CoV-2.

Improved Neutralization of Omicron BA.4/5, BA.4.6, BA.2.75.2, BQ.1.1, and XBB.1 with Bivalent BA.4/5 Vaccine (preprint)

The BNT162b2 bivalent BA.4/5 COVID-19 vaccine has been authorized to mitigate COVID-19 due to current Omicron and potentially future variants. New sublineages of SARS-CoV-2 Omicron continue to emerge and have acquired additional mutations, particularly in the spike protein, that may lead to improved viral fitness and immune evasion. The present study characterized neutralization activities against new Omicron sublineages BA.4.6, BA.2.75.2, BQ.1.1, and XBB.1 after a 4th dose (following three doses of BNT162b2) of either the original monovalent BNT162b2 or the bivalent BA.4/5 booster in individuals >55 years of age. For all participants, the 4th dose of monovalent BNT162b2 vaccine induced a 3.0, 2.9, 2.3, 2.1, 1.8, and 1.5 geometric mean neutralizing titer fold rise (GMFR) against USA/WA1-2020 (a strain isolated in January 2020), BA.4/5, BA.4.6, BA.2.75.2, BQ.1.1, and XBB.1, respectively; the bivalent vaccine induced 5.8, 13.0, 11.1, 6.7, 8.7, and 4.8 GMFRs. For individuals without SARS-CoV-2 infection history, BNT162b2 monovalent induced 4.4, 3.0, 2.5, 2.0, 1.5, and 1.3 GMFRs, respectively; the bivalent vaccine induced 9.9, 26.4, 22.2, 8.4, 12.6, and 4.7 GMFRs. These data suggest the bivalent BA.4/5 vaccine is more immunogenic than the original BNT162b2 monovalent vaccine against circulating Omicron sublineages, including BQ.1.1 that is becoming prevalent globally.

Cross-neutralization and viral fitness of SARS-CoV-2 Omicron sublineages (preprint)

The rapid evolution of SARS-CoV-2 Omicron sublineages mandates a better un-derstanding of viral replication and cross-neutralization among these sublineages. Here we used K18-hACE2 mice and primary human airway cultures to examine the viral fit-ness and antigenic relationship among Omicron sublineages. In both K18-hACE2 mice and human airway cultures, Omicron sublineages exhibited a replication order of BA.5 [≥] BA.2 [≥] BA.2.12.1 > BA.1; no difference in body weight loss was observed among differ-ent sublineage-infected mice. The BA.1-, BA.2-, BA.2.12.1-, and BA.5-infected mice de-veloped distinguisable cross-neutralizations against Omicron sublineages, but exhibited little neutralization against the index virus (i.e., USA-WA1/2020) or the Delta variant. Surprisingly, the BA.5-infected mice developed higher neutralization activity against heterologous BA.2 and BA.2.12.1 than that against homologous BA.5; serum neutralizing titers did not always correlate with viral replication levels in infected animals. Our results reveal distinct antigenic cartography of Omicron sublineages and support the bivalent vaccine approach.

Low neutralization of SARS-CoV-2 Omicron BA.2.75.2, BQ.1.1, and XBB.1 by 4 doses of parental mRNA vaccine or a BA.5-bivalent booster (preprint)

The newly emerged SARS-CoV-2 Omicron BQ.1.1, XBB.1, and other sublineages have accumulated additional spike mutations that may affect vaccine effectiveness. Here we report neutralizing activities of three human serum panels collected from individuals 1-3 months after dose 4 of parental mRNA vaccine (post-dose-4), 1 month after a BA.5-bivalent-booster (BA.5-bivalent-booster), or 1 month after a BA.5-bivalent-booster with previous SARS-CoV-2 infection (BA.5-bivalent-booster-infection). Post-dose-4 sera neutralized USA-WA1/2020, BA.5, BF.7, BA.4.6, BA.2.75.2, BQ.1.1, and XBB.1 SARS-CoV-2 with geometric mean titers (GMTs) of 1533, 95, 69, 62, 26, 22, and 15, respectively; BA.5-bivalent-booster sera improved the GMTs to 3620, 298, 305, 183, 98, 73, and 35; BA.5-bivalent-booster-infection sera further increased the GMTs to 5776, 1558,1223, 744, 367, 267, and 103. Thus, although BA.5-bivalent-booster elicits better neutralization than parental vaccine, it does not produce robust neutralization against the newly emerged Omicron BA.2.75.2, BQ.1.1, and XBB.1. Previous infection enhances the magnitude and breadth of BA.5-bivalent-booster-elicited neutralization.

Determinants and Mechanisms of the Low Fusogenicity and Endosomal Entry of Omicron Subvariants (preprint)

The rapid spread and strong immune evasion of the SARS-CoV-2 Omicron subvariants has raised serious concerns for the global COVID-19 pandemic. These new variants exhibit reduced fusogenicity and increased endosomal entry pathway utilization compared to the ancestral D614G variant, the underlying mechanisms of which remain elusive. Here we show that the C-terminal S1 mutations of the BA.1.1 subvariant, H655Y and T547K, critically govern the low fusogenicity of Omicron. Notably, H655Y also dictates the enhanced endosome entry pathway utilization. Mechanistically, T547K and H655Y likely stabilize the spike trimer conformation, as shown by increased molecular interactions in structural modeling as well as reduced S1 shedding. Importantly, the H655Y mutation also determines the low fusogenicity and high dependence on the endosomal entry pathway of other Omicron subvariants, including BA.2, BA.2.12.1, BA.4/5 and BA.2.75. These results uncover mechanisms governing Omicron subvariant entry and provide insights into altered Omicron tissue tropism and pathogenesis.

Development of highly potent non-covalent inhibitors of SARS-CoV-2 3CLpro (preprint)

The SARS-CoV-2 virus is the causal agent of the ongoing pandemic of coronavirus disease 2019 (COVID-19). There is an urgent need for potent, specific antiviral compounds against SARS-CoV-2. The 3C-like protease (3CLpro) is an essential enzyme for the replication of SARS-CoV-2 and other coronaviruses, and thus is a target for coronavirus drug discovery. Nearly all inhibitors of coronavirus 3CLpro reported so far are covalent inhibitors. Here, we report the development of specific, non-covalent inhibitors of 3CLpro. The most potent one, WU-04, effectively blocks SARS-CoV-2 replications in human cells with EC 50 values in the 10-nM range. WU-04 also inhibits the 3CLpro of SARS-CoV and MERS-CoV with high potency, indicating that it is a pan-inhibitor of coronavirus 3CLpro. WU-04 showed anti-SARS-CoV-2 activity similar to that of PF-07321332 (Nirmatrelvir) in K18-hACE2 mice when the same dose was administered orally. Thus, WU-04 is a promising drug candidate for coronavirus treatment. One-Sentence Summary A oral non-covalent inhibitor of 3C-like protease effectively inhibits SARS-CoV-2 replication.

Integrated multi-omics analyses identify key anti-viral host factors and pathways controlling SARS-CoV-2 infection (preprint)

Host anti-viral factors are essential for controlling SARS-CoV-2 infection but remain largely unknown due to the biases of previous large-scale studies toward pro-viral host factors. To fill in this knowledge gap, we performed a genome-wide CRISPR dropout screen and integrated analyses of the multi-omics data of the CRISPR screen, genome-wide association studies, single-cell RNA-seq, and host-virus proteins or protein/RNA interactome. This study has uncovered many host factors that were missed by previous studies, including the components of V-ATPases, ESCRT, and N-glycosylation pathways that modulated viral entry and/or replication. The cohesin complex was also identified as a novel anti-viral pathway, suggesting an important role of three-dimensional chromatin organization in mediating host-viral interaction. Furthermore, we discovered an anti-viral regulator KLF5, a transcriptional factor involved in sphingolipid metabolism, which was up-regulated and harbored genetic variations linked to the COVID-19 patients with severe symptoms. Our results provide a resource for understanding the host anti-viral network during SARS-CoV-2 infection and may help develop new countermeasure strategies.

Neutralization of SARS-CoV-2 Omicron sublineages by 4 doses of mRNA vaccine (preprint)

Since the initial emergence of SARS-CoV-2 Omicron BA.1, several Omicron sublineages have emerged, leading to BA.5 as the current dominant sublineage. Here we report the neutralization of different Omicron sublineages by human sera collected from individuals who had distinct mRNA vaccination and/or BA.1 infection. Four-dose-vaccine sera neutralize the original USA-WA1/2020, Omicron BA.1, BA.2, BA.212.1, BA.3, and BA.4/5 viruses with geometric mean titers (GMTs) of 1554, 357, 236, 236, 165, and 95, respectively; 2-dose-vaccine-plus-BA.1-infection sera exhibit GMTs of 2114, 1705, 730, 961, 813, and 274, respectively; and 3-dose-vaccine-plus-BA.1-infection sera show GMTs of 2962, 2038, 983, 1190, 1019, and 297, respectively. Thus, 4-dose-vaccine elicits the lowest neutralization against BA.5; 2-dose-vaccine-plus-BA.1-infection elicits significantly higher GMTs against Omicron sublineages than 4-dose-vaccine, and 3-dose-vaccine-plus-BA.1-infection elicits slightly higher GMTs (statistically insignificant) than the 2-dose-vaccine-plus-BA.1-infection. Our results support the inclusion of the BA.5 spike for future vaccine booster doses.

Neutralization of Omicron sublineages and Deltacron SARS-CoV-2 by 3 doses of BNT162b2 vaccine or BA.1 infection (preprint)

Distinct SARS-CoV-2 Omicron sublineages have evolved showing increased fitness and immune evasion than the original Omicron variant BA.1. Here we report the neutralization activity of sera from BNT162b2 vaccinated individuals or unimmunized Omicron BA.1-infected individuals against Omicron sublineages and Deltacron variant (XD). BNT162b2 post-dose 3 immune sera neutralized USA-WA1/2020, Omicron BA.1-, BA.2-, BA.2.12.1-, BA.3-, BA.4/5-, and XD-spike SARS-CoV-2s with geometric mean titers (GMTs) of 1335, 393, 298, 315, 216, 103, and 301, respectively; thus, BA.4/5 SARS-CoV-2 spike variant showed the highest propensity to evade vaccine neutralization compared to the original Omicron variants BA.1. BA.1-convalescent sera neutralized USA-WA1/2020, BA.1-, BA.2-, BA.2.12.1-, BA.3-, BA.4/5-, and Deltacron-spike SARS-CoV-2s with GMTs of 15, 430, 110, 109, 102, 25, and 284, respectively. The low neutralization titers of vaccinated sera or convalescent sera from BA.1 infected individuals against the emerging and rapidly spreading Omicron BA.4/5 variants provide important results for consideration in the selection of an updated vaccine in the current Omicron wave.

Antibody escape and cryptic cross-domain stabilization in the SARS CoV-2 Omicron spike protein (preprint)

The worldwide spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has led to the repeated emergence of variants of concern. The Omicron variant has two dominant sub-lineages, BA.1 and BA.2, each with unprecedented numbers of nonsynonymous and indel spike protein mutations: 33 and 29, respectively. Some of these mutations individually increase transmissibility and enhance immune evasion, but their interactions within the Omicron mutational background is unknown. We characterize the molecular effects of all Omicron spike mutations on expression, human ACE2 receptor affinity, and neutralizing antibody recognition. We show that key mutations enable escape from neutralizing antibodies at a variety of epitopes. Stabilizing mutations in the N-terminal and S2 domains of the spike protein compensate for destabilizing mutations in the receptor binding domain, thereby enabling the record number of mutations in Omicron sub-lineages. Taken together, our results provide a comprehensive account of the mutational effects in the Omicron spike protein and illuminate previously unknown mechanisms of how the N-terminal domain can compensate for destabilizing mutations within the more evolutionarily constrained RBD.

SARS-CoV-2 Delta breakthrough infections in vaccinated patients (preprint)

The continuous emergence of SARS-CoV-2 variants with increased transmission and immune evasion has caused breakthrough infections in vaccinated population. It is important to determine the threshold of neutralizing antibody titers that permit breakthrough infections. Here we tested the neutralization titers of vaccinated patients who contracted Delta variant. All 75 patients with Delta breakthrough infections exhibited neutralization titers (NT50) of less than 70. Among the breakthrough patients, 76%, 18.7%, and 5.3% of them had the NT50 ranges of <20, 20-50, and 50-69, respectively. These clinical laboratory results have implications in vaccine strategy and public health policy.

Cross neutralization of Omicron BA.1 against BA.2 and BA.3 SARS-CoV-2 (preprint)

The Omicron SARS-CoV-2 has three distinct sublineages, among which sublineage BA.1 is responsible for the initial Omicron surge and is now being replaced by BA.2 worldwide, whereas BA.3 is currently at a low frequency. The ongoing BA.1-to-BA.2 replacement underscores the importance to understand the cross-neutralization among the three Omicron sublineages. Here we tested the neutralization of BA.1-infected human sera against BA.2, BA.3, and USA/WA1-2020 (a strain isolated in late January 2020). The BA.1-infected sera neutralized BA.1, BA.2, BA.3, and USA/WA1-2020 SARS-CoV-2s with geometric mean titers (GMTs) of 445, 107, 102, and 16, respectively. Thus, the neutralizing GMTs against heterologous BA.2, BA.3, and USA/WA1-2020 were 4.2-, 4.4-, and 28.4-fold lower than the GMT against homologous BA.1, respectively. These findings have implications for vaccine strategy.

Neutralization of Omicron BA.1, BA.2, and BA.3 SARS-CoV-2 by 3 doses of BNT162b2 vaccine (preprint)

The newly emerged Omicron SARS-CoV-2 has 3 distinct sublineages: BA.1, BA.2, and BA.3. BA.1 accounts for the initial surge and is being replaced by BA.2, whereas BA.3 is at a low prevalence at this time. Here we report the neutralization of BNT162b2-vaccinated sera (collected at 1 month after dose 3) against the three Omicron sublineages. To facilitate the neutralization testing, we engineered the complete BA.1, BA.2, or BA.3 spike into an mNeonGreen USA-WA1/2020 SRAS-CoV-2. All BNT162b2-vaccinated sera neutralized USA-WA1/2020, BA.1-, BA.2-, and BA.3-spike SARS-CoV-2s with titers of >20; the neutralization geometric mean titers (GMTs) against the four viruses were 1211, 336, 300, and 190, respectively. Thus, the BA.1-, BA.2-, and BA.3-spike SARS-CoV-2s were 3.6-, 4.0-, and 6.4-fold less efficiently neutralized than the USA-WA1/2020, respectively. Our data have implications in vaccine strategy and understanding the biology of Omicron sublineages.